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probability of K symbol errors in a block of n m-bit symbols. 
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SECTION I 


THE MODEL 

In estimating the performance of an error-correcting device on 
a specific communication channel, it is necessary to find a meaning¬ 
ful, yet tractable, mathematical model for that channel. An examination 
of data from real channels suggests that most channels pass through 
three distinct phases. The first phase, which nearly any error- 
correcting scheme can handle successfully, is that of long periods of 
practically error-free transmission. The second phase, which is the 
antithesis of the first phase in that no scheme can expect to correct 
it, is that of complete loss of signal for substantial periods of time. 
The third phase might be described as a generalized sputter or bursts 
of error bursts, and it is this phase, if it occurs frequently, for 
which error correctors should be designed. This third phase may 
sometimes be described by means of a two-state model. The three 
phases suggest a Markov process with four states, but such a process 
is mathematically unwiedly. However, when the three phase picture is 
reasonably correct and the sputter phase is accurately modeled for 
purposes of estimating coder performance, it is a simple matter to 
make corrections in such estimates for those periods when transmission 
is nearly error free or when the signal is lost. 

The use of such a two-state model was suggested to us by some work 
by Gilbert [1]. Gilbert describes a model for binary error distributions 
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in channels subject to noise bursts. Let {x.} be the error process with 

x # = 1 for an error in the i^ demodulated bit, x. = 0 for no error, 
i l 

Two states, designated G and B, of the channel are postulated such that 

th 

at the i bit the state S. is G or B, and the state S 4| , at bit i + 1 

l l+l 

depends only on Thus, the state sequence {S^( is a simple Markov 

chain described by the two transition probabilities, P: G -> B and 
p B -» G. We also use Gilbert's notation Q = 1 - P and q = 1 - p for 
the G -* G and B -* B transitions respectively Let h,k denote respectively 
the probabilities of a correct bit in B and in G Then P[x^=l( = 1 - h - h 1 
if S. - B and P{x.=l] - 1 - k - k* if S. - G 

l l i 

Using this model with k = 1, Gilbert obtained a good fit to certain 

phone line error data. The statistic fitted was the probability of 

occurrence of zero (= error free) runs of length at least K (The use 

of k < 1 has also been considered by Elliott [5].) 

At MITRE, we have found sets of values for the model parameters 

P, p, h, k, which yield good fits to error data for several different 

types of communication media One important statistic involved is the 

probability of specific error densities in various length blocks. 

Given that we have p, P, h,k, [2,3,4], we present the recursive 

method used to calculate the probability of K symbol errors in a block 

of n symbols where each symbol consists of m bits. 
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SECTION II 


RECURSIONS - PART I 

In Appendix II of [2] (or Appendix B of [3] or [4]) we presented a 
brief outline of the recursive method used to calculate the probability 
of K symbol errors in a block of n symbols where each symbol consists 
of m bits Now we present the outline in full. Since the last 
documentation, the technique has been extended to permit the n-symbol 
blocks to be interleaved or time spread. (To interleave s blocks means 
to transmit sequentially the first symbol of each of s blocks followed by 
the second symbol of each of those blocks, etc. Thus, a given m-bit 
symbol is transmitted s symbol times after the symbol preceding it in 
its block.) 

Letting T and U represent either of the states G and B, we define 
TOU(t) = P(x x ■ • • • ■ ■ 0 and S^_ = U|S Q = T ) 

TlU(t) = P(for some i £ t, x ■ 1 and c U|Sq = T) 

Then 


G0G(1) = Qk 

G0B(1) = Ph 

G1G(1) = Qk' 

G1B(1) = Ph' 

B0B(1) = qh 

B0G(1) = pk 

B1B(1) = qh' 

B1G(1) = pk' 


and 

GOG(t) = [GOB(t-1) p + GOG(t-l) Q] k 

GlG(t) = [GOB(t-l) p + GOG(t-l) Q] k' + GlB(t-l) p + GlG(t-l) Q 
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GOB(C) = [GOB(t-l) q + GOG(t-l) P] h 

GlB(t) = [GOB(t-1) q + GOG(t-l) P] h' + GlB(t-i; q + GlG(t-l) P 

BOB(t) = [BOB(t-l) q + BOG(t-l) P] h 

BlB(t) = [BOB(t-l) q + BOG(t-l) P] h' + BlB(t-l) q + BlG(t-l) P 

BOG(t) = [BOG(t-l) Q + BOB(t-l) p] k 

B1G(C) = [BOG(t-l) Q + BOB(t-l) p] k' + BlB(t-l) p + BlG(t-l) Q 

Shortly after the program was written, we discovered that GOG(m), 
GOB(m), GlG(m), GlB(m), BOG(m), BOB(m), BlG(m), BlB(m) could be obtained 
from a difference equation in powers of J and L (see below). Since on 
our computer (IBM 7030), it took under a second to compute all eight 
quantities, we decided not to use the difference equation. However, we 
work out two, and give all eight results. 

GOG(t) = [GOG(t-l) Q + GOB(t-1) p} k 
GOB(t) = [GOB(t-l) q +G0G(t-l) P} h 

The eigenvalues come from the 2nd order linear difference equation: 

f t+l ‘ (Qk + qh) f t ' (P ' Q) f t-l = °* 
that is, 2j = Qk + qh W(Qk + qh) 2 + 4hk(p-Q) 

and 21. = Qk + qh - ^(Qk + qh)^ + 4hk(p-Q) 

Thus we have 

GOG(t) - ^ + 0^ L fc , 

GOB(t) = fi l l' 1 
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and we get 


<*!> 0 ' 2 , ^ and ^ 


from the initial conditions 

GOG(O) = 1, G0G(1) = Qk 
and GOB(O) = 0, G0B(1) = Ph . 

So a l + °2 = + O^L = Qk, which yields 

GOG(t) = {(Qk - L)/(J - L)} J l + {(J - Qk)/(J - L)} L* 

Also, + ^2 = ^1 J + ^2 L = Ph> which y ields 

GOB(t) = {Ph/(J - L)} (J C - L C ) 

All eight solutions are as follows: 

GOG(m) = {(Qk - L)/(J - L)} J™ + {(J - Qk)/(J - L)} L™ 

GOB(tn) = {Ph/(J - L)} (J™ - L™) 

BOG(m) = [pk/(J - L)} (J™ - L™) 

BOB(tn) = {(qh - L)/(J - L)} J™ + {(J - qh)/(J - L)} L™ 

GlG(tn) = p/(p + P) + [P/(p + P)} (Q~p) m - {(J - qh)/(J - L)} J 1 " 
- {(qh - L)/(J - L)} L m 

GlB(tn) - ^ [1 - (Q~P) m ] - flI ^ 

BlB(m) = —j + (Q-p) m - ^ + !Qk-LjL m ] 

BlG(tn) = [1 - (Q-P)™] - ^ (J 1 " - L™) 
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SECTION III 


RECURSIONS - PART 2 


Again letting T and U represent either of the states G and B, 


we define 



x 


- 0,S 


sm 



T) 


= P(m-bit symbol after s symbol times is correct and 

ends in state ujstate T) 


TlUI(s) = P(for some 1 £ i £ m, x / ...=1,S = u|S n = T) 

(s-I)mfi sm 1 0 


= P(m-bit symbol after s symbol times has at least one bit 

error and ends in state u|state T) 


Let GXG = GOG(m) + GlG(m), GXB - GOB(m) + GlB(m), BXG = BOG(m) + BlG(m), 
and BXB = BOB(m) + BlB(m) 

Then GOGI(s) = GXG GOGI(s-l) + GXB * BOGI(s-l) and BOGI(s) = 

BXG GOGI(s-l) + BXB * BOGI(s-l). 

(There will be similar equations in G1BI, B1BI and GOBI, BOBI, and 

G1BI, B1B1, but they have the same eigenvalues and will be omitted.) 

TOGI(s) - [GXG + BXB] TOGI(s-l) + [GXG*BXB-GXB BXG] T0GI(s-2) « 0 

This yields the eigenvalues: 

2 CT = GXG + BXB + >/[GXG + BXB] 2 -4 [GXG • BXB - GXB ■ BXG] 

2 t = GXG + BXB -+ BXB] 2 -4 [GXG ' BXB - GXB • BXG] 



6 




Since GXG + GXB = BXB + BXG = 1, 


a T = GXG • BXB - GXB * BXG = (1-GXB) (1-BXG) - GXB • BXG = 1-GXB-BXG 

a + T = GXG + BXB = 2 - GXB - BXG = 1 + a T . 

Hence a = 1, and T = 1 - GXB - BXG = GXG + BXB - 1. 

Thus we have 

GOGI(s) = X 1 1 S + \ 2 t S 
and BOGI(s) = 1 S + 

and we get and from the initial conditions. 

G0GI(1) = GOG(m), G0GI(2) = GXG * GOG(m) + GXB • BOG(m) 

BOGI(l) = BOG(m), B0GI(2) = BXG * GOG(m) + BXB • BOG(m) 

Hence 

A 1 + X 2 ' T = G0G ( m ^> \ ’ l 2 + ’ T 2 = GXG * G0G ( m ) + GXB * BOG(m). 

Solving, X = [GOG(m) (1 - BXB) + BOG(m) • GXB] / (1 - t) 

\ 2 = [GOG(m) BXB - BOG(m) GXG] / (1 - t) = [GOG(m) - ] / t 

Also, 

H • 1 + ^2 ‘ T = BOG(tn), ^ • l 2 + ^ • T 2 = BXG * GOG(m) + BXB • BOG(m). 

Solving, = [BOG(m) (1 - GXG) + BXG * GOG(m)] / (1 - t) 

42 = [BOG(m) - / T 
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Therefore, GOGI(s) = X^ ' I s + X^ * t S , and BOGI(s) = Mj'l S + ti 2 * T S , 
where X^, are given above. Since s is fixed for any given 

application, we will refer to GOGI(s) as GOGI, and to BOGl(s) as BOGI. 

Consider the first i m-bit symbols of a random interleaved block 
Let GB(i,j) = P(j symbol errors in i symbols and b|Sq * G). 

Similarly, we define GG(i,j), BB(i,j) and BG(i,j). 

Then 


GG(1,0) - GOGI 

GB(1,0) - GOBI 

GG(1,1) = G1GI 

GB(1,1) = G1BI 

BG(1,0) = BOGI 

BB(1,0) = BOBI 

BG(1,1) = BlGI 

BB (1,1) = BlBI 


Finally, for i = 2, ■ , n and j =0, ' ' , i 

GG ( i , j ) = GG ( i -1 , j .) GG (1,0) + GB(i-l,j) BG(1,0; 

+ GG(i-l.j-l) GG(1,1) + GB(i-1,j-1) BG(1,1) 

GB(i , j) = GG(i-1,j) GB(1,0) + GB(i-l,j) BB(1,0J 

+ GG(i-1,j-1) GB(1,1) + GB(i-1,j-1) BB(1,1) 


BG(i , j) = BG(i-l,j) GG(1,0) + BB(i-l,j) BG(1,0) 

+ BG(i-l,j-1) GG(1,1) + BB(i-1,j-1) BG(1,1) 

BB(i , j) = BG(i-1,j) GB(1,0) + BB(i-l,j) BB(1,0) 

+ BG(i-l,j-1) GB(1 ,1) + BB(i-l , j) BB(1,1). 
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Finally, 


P(random bit is in G) = a. = and hence 


P(K symbol errors in n symbols with s blocks interleaved) 


CL [GG(n,K) + GB(n,K)] + (1 - a) [BG(n,K) + BB(n,K)] 


9 


SECTION IV 


INPUT FOR THE PROGRAM 


XM 

XN 

XNEST 


XIPER 

IK 

CP 

SP 

H 

SK 


No of bits/symbol 
No. of symbols/block 

Largest number of symbol errors to be considered 
(if this field is blank, XNEST = XN) 

Number of interleaved symbols (0 or 1 means 1) 

0 or not equal to 0 (0 means continue with CP, SP, H, SK; 
not equal to 0 means read new parameters) 

P 

P 

h 

k 


» 


« 
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SECTION V 


OUTPUT OF THE PROGRAM 

TIMEX = A8 representation of the time read from IBM 7030 Time Clock 
by a STRAP coded routine (one can call his routine or omit 
it altogether) 

CP, SP, H, SK as above 

ALPHA = P (random bit is in G) 

M (or XM) as above 

N (or XN) as above 

NS = No. of recursion terms to be attempted 
WPI = P (symbol error) 

WY155 = P (no errors when s = 1) 

WPMU2 = mean number of errors in a block 
FCMEAN = mean number of errors given an error occurred (when s = 1) 
For J, L, A, B see [2,3,^] 

IS = number of errors 
P * P (IS symbol errors in a block) 

R = P (IS symbol errors in a block|error occurred) 

Q * P (£ IS symbol errors in a block) 

QH * P (> IS symbol errors in a block) 

S * mean number of bits between blocks with > IS symbol errors 

PBAR * contribution to mean number of errors per block made by 

probabilities actually calculated (if one wants the whole 
mean, then NEST = 0 or N; the same applies to VAR, SVAR, 
and CMEAN) 
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VAR 

SVAR 

CMEAN 


contribution to variance about PBAR made by probabilities 
actually calculated by the recursion 

= approximate standard deviation 

= contribution to mean number of errors given an error 
occurred in the block made by probabilities actually 
calculated by recursion 
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APPENDIX 


PROGRAM 


V fcftSI dH 55 72R /G5 - CORRECTION LE¥Etl £T3/2B?i5 


00000 

0.QU0I 


C...T0 CALCULATE the recursion probabilities 
COMMUN 771/ TIMEX 

DIMENSION GGI515)t G6{515), BGj515_L» OBC515) 





00005 

00006 


_ _ __ ■ -^S fl l2Crt l , G1BI2031 = 

X H l50I,Wl5ai, XNESI (5 0 1* n=T^ XTF*iU 50) 

404 REA01701, LK 

1701 FORMAT ill) 

C...KEAO IN THE NO. OF BITS PER SYMBOL, THE NO. OF SYMBOLS, 




jswmr~ 7gr 


00009 

00010 . 


C.. .R 

604 

702 .. 


HO* MANY SYMBOLS TO GO THRU * HOk MANY 
Ll-tr 

T EQgMAT I IMX, 16, 2 X t 13, . _ . .141= 

EAD IN CP, SP, H, SK 
READ 702, IK, CP, SP, H, SK 
FORMAT (11, E17.10, 3CE18.il)) 


iTitKi 


m^z 


m Z&& K*l, LK 


M XMTK T 

| XNIK) 


00014 

00016 

00016 


NEST - XNEST(K) 
IPERD =X IPER(K) 
IF(NEST.EQ.O) NEST 


000 IT 
0B31R 


,DQSTS 


NS * NEST + 1 

CALL TIME 

PRINT 300, TIMEX 


00020 


00021 


300 


00025 

00026 

00027 


FORMAT(1H1,Ad ,21X, 77H*CALCULAT1UN OF GILBERT CHANNEL ERROR PROBAB 
llLlTItS USING RECURSION FORMULAS*/////) 

MUD _=_M*N______ 

fMUO - FLOAT (MUD) 

MU 01 HUD - 1 

_=^ !> - H _ 1_ . _ . __ 

SKCUMP = 1* - SK 

CQ = 1 * - CP 


so 


1 • - SP 


GAMMA = 1 ./ 1 CP 4 -SP) 
ALPHA » SP*GAMMA 
BETA = CP*GAMMA 


00031 

00032 

00033 


301 


PI = (SP* SKCUMP ♦ CP*HCOMP)*GAMMA 
PRINT 301, CP, SP, H, SK, ALPHA 
FORMAT (5X» 2 7H0HANNEL PARAMETERS 


CP = ,E11•6, 5X, 5HSP * 


4ElU64X f H * f F8.6, 5X, 5H K = ,Elfc*lG*S3U mALPHA = * FBJA) 

PRINT 302* Mt N* PiSi Pi* lPEfUl 

ZfitfflE F AR&fgT SRS- = M-= +-ZW+ -1 31IT W " 


112X, 35HN0. OF RECURSION TERMS REQUESTED = , 

2 Eli.7, 75X, 14, 

3 $,/) 


15/ 7X, SP1 = S, 

$ BLOCKS INTERLEAVED 
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* 


0CC36 

00037 

00038 

mN 

WT 

HR 

* 

SP - CQ 

8K*CQ + H* SO 

WT**2 € 4.*H*SK*WN 


f , 

00039 

00040 

MSUKT 
w J = 

= SORT(HR) 

t *T ♦ HSOR T)/2. 



00041 

00U42 

WL 

WA 


(WT - wS0RT)/2. 

WJ ♦ wN *(CP*Sk*HCURP + SP*H*SKCOMP ) / 




1 (CP*HCUMP + *P*$KCUMP ) 


UQC43 WA = HA / WSQRT 


00044 

W 8 = 

HA - 1 . 

00u45 

H X 2 = 

(WA*(1* —WJ**M J/(1. —W J J) - (w B* ( l.-WL**M)/(l.-WL) ) 

0 0 0 4 6 

00G4 7 
00048 

HP I = 
HPMU2 

K J * 

PI * HX2 
= FLOAT(N) * HP1 

F LGAT t MUO) * ALGGlO(WJ) 


0UG49 KL - FLUAT(MUD) * ALUGlOIwU 


00050 

1FIR0. 

LT. 

-50. ) 

F 

= 0* 

00051 

IFIRJ. 

GE. 

-50.) 

F 

= WJ**MUO 

000 52 

( F (R L . 

LT. 

-60. ) 

G 

= 0. 

00053 

(Ft*L . 

GE. 

-50. ) 

G 

= WL**MUL> 


00054_ HY155 = 1. -( trtA* { 1 .-F ) / t 1 .-HJii-HdM 1 .-G) / 1 U-HL ) ) *P 1 

00055 fCrtfcAW = WPMU2 / (1. - WY155) 

00056 PRffSiT 806, WPJ, W Y 1 5 5 » UPHU2 ,FbMfc AN, W J , *L* WA, Wtt 

00057_ 806 FORMAT (IX, 20HPRQ6 SYMBOL ERROR = ,£12,6* 1H*, 2X* I7HPRQB jNO £R 

1K0RS = , b 1 2.6 * 1H*, 2 X, OH Mt AN = ,tl2.6, 1H*, 2X, SCUNOITlUNAL ME 

2AN = i, 11 2 • 6 , $*$ / 1JX, 4HJ = , 3 X, £16.10, 1H*, 

31 5 X, <♦ HL = , 3X, F9.6, 1H* , 6X , 4H A > , 3X, F9 .6, lH*, 20X, 4HB « 

4fF9.6t 1H*///1 


00058 

GC059 

G 0 C 6 C 

00061 

00062 

00063 

00064 

00065 

00066 

00067 

00068 

0006* 
00070 
00071 
00072 
00073 
00074 
£00075 
00076 
00077 
00078 
0 00 79 
00080 
00081 


GOG(1) 
GlG(l) 
GOB(1) 
Glbl11 
BJBl 1J 
BIB(1> 
BOG(1) 
B1G(1) 

I Ft M.EU 
00 50 
J1 = J- 


= CO* SK 
* €0 - 
= C P* H 
= CP - 
= SQ*H 
= SQ - 
= SP* SK, 

= SP - 
li GO TU 
»M 


G0G(1> 


GOB(1J 
BOB(1) 

BOG(1) 

8u 


50 

80 


TA 

GOGtOJ * 
GIG(J) = 
TB 

GOBtJ) = 
G1B (J) « 
T£ 

BOB(J) = 
B1BU1 = 
TO 

BOGtJ) = 
B1G(J) = 
PRINT 3C3 


* G1G(JU*CQ 
G OG t J1) * CP 


* Glri(Jil*5P 


J = 2 , 

1_ ...___ 

= G0b("J1)*SP + G0GIJ1)*C0 
» T A* SK 
= TA - GOG(J) 

= GOBt J1)*S0 < 

= T B*H 

= TB - GOB(J) 

= BQd(Jl)*SO ♦ BOG(Jl)*CP 
= TC*H 

TC - BOB(J) 

BOG(Jl>*CW 
TO* SK 

TO - BOG(J) ♦ BIG ( J1 > *C0 + BlbtJi ) *S P 
GOG(M), GlG(H), GOB(Ml * GlB(rt) 


♦ GIBtJl)*SQ + G1G(Jl)*CP 


* BIB(31)*S0 
B OB( JlJ * SP 
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00092 303 fOKMAt 6HG0G ? , bl2«5* 5X? 6Hult» = * £I2»5, 5X^ 6MG08 — y - 

=-^=-=^ SXg^ttOig_-i«> ♦ €12,3* _ - -- ~- ^== 

000B3 PRINT 304, BOG(M) , BIG(M), B08(M), BIB(M) 

000B4 304 FORMAT (5X, 6HB0G = , E12.5, 5X f 6HB1G = , E12.5, 5X, 6HB0B * * 

IE 12•5, 5X f 6HBIB = , E12.5//I 

00085 ^^^ ^TR PI1 - ALPHA* (GlSTMf * G1G1H) ) * dETA»( BIB i *i + 0101*0 

00086 . PRINT 505, RPII 

QOOB7 — — 505 fO^HAf JiX* 30HPROB RANDOM SVM80L IN ERROR ■ , EL2«6, 1M+///} 

OOOBB CALL TIME 

000B9 GG(1) = GOG(M) 

00090 GG(2) = G1G(M) 

80091 GBC1I - G08f«l=^=7== 

jSg6« — g*t2> - 

00094 BG(2) = BIG(M) 

00095 BB< 1) = BOB ( M) 

00096 BB(2) * B1B(M) 

00092- IF (fPEKD ,N£* 01 GO TO 9001 

00098 GOG I = bOG f MI 

00099 GlCI * GIGIMI 

00100 GOBI = GOB(M) 

00101 GIB I = GIB(M) 

00102 _BOG I = BOG(M) __ _ 

00103 BIG! = BiGIM) 




00107 9001 GXG = bGU) ♦ GGt 21 
00108 GXB ^ GB(I> 4 Gb(2 I 
00109 _ 6XG = 6GI1 ) * BG I 21 



00119 PRINT 9003, SIG, TAU 
00120 9003 FORMA! (6X, SSIGMA = 
-00121 AG* ( GG< 1) *GT* 


$, E20.12, 10X, STAU 
GXB * BG(1)) / DEN 


$, E20.12/' 


mx ar 

-«*» 

80124 

00125 

00126 

00127 


GOGI = AG 
BOGI = AB 
AG = (GG(2) 


*GT 


CG 

CB 

♦ 


♦ G 

♦ G 
GXB 


♦ BG ( 2)) / DEN 


00129 


AB 


fGGf21 - AG I / TAU 

10X0 * 00£2t * BGI2I * 8T) / 0£H 


% 
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0G130 

00131 

00132 

0 C 1 33 

0013* 

00135 

00136_ 

0C137 

0013B 

00139 

001*0 

0C141 

00142 

001*3 

00144 

001*5 


CB = IBg(2) - AB) / TAU 
G1GI = Ab ♦ CG * G 

ttlGI = AB ♦ CB * G 

AG_= <_GB11) *GT ♦ GXb 

.CG = (GBill-AG) / TAU 
AB = (BXC » Will + BttUl 
CB = IBBi1) - AB) / TAU 


* B6<1)1 / OEN 


* BT) / O £H 


OEN 


CG * G 
C B * G 

*6T _♦_GXB ♦ BB(2 )) / 

CG = (G8( 2) - AG) / TAU 

AB * <8XG * Gbl2) ♦ Bb(2) * BU / D£N 

CB = IBB12) - AB) / TAU 


GCBI = Ab 
BOBI = AB 
AG = (gB(2) 


00146 


00161 


G1BI = Ag ♦ Co * G 

B1BI = AB ♦ CB * G 

_ PRINT 9003 t GOG1t G1G1> GOBI t GlBlt BOG 1» B 1GI, BOB It B1B1 _ 

900!* £GitMAT I/sX* 7MG0GI = W E12^5* 4X* 7HG1GI , £12.5* 4X, 7HGU&I 

i €12.5« *X» 7HG10I * t £12.5/ 5X f 7HB0G1 • , £12.5* *X f 7 hB1gI 

7H30BI = , £12.5, *X, 7HBIS1 = , E12.5//1 


Z £12.5* 4Xi 


p -*! 


00147 

C 

00146 

9002 1F(NEST.EQ. 

.. .CALCULATE THE 
DO 51 L1=2 

00149 

LZ = 

MING 

00150 

LY = 

Cl 

00151 

IFIL1 

.GT. 1 

00152 

IFILI 

.6T. 1 

00153 

N1 = 

LI ♦ 1 

00154 

GG(N1) 


» N£$T + 1 


Go IL 1)*61GI 


GB(LI)*B1G( 


00155 G&IN1) 

00156 BGiNU 

00157 BBlril) • 

00158 97 CONTINUk 

00159 DU 02 L2=2,LY 

00160 N2 = LY - L2 


GG(LI) *61BI 
BG < L11*u1 gI 
BGIL1)*6 IB I 


GSTL1)*B1BI 

MTi L U^GI 

BBt€n*BlBl 


+ 2 


N3 = N2 - 1 


00162 YAH1 *GGIN2)*G0Gi 

00163 GB4N2)=6BINZ)*BOBI 

00164 G&TN2) = YAH1 

00165 YAH2 =BG(N2)*G0GI 

00166 BB(N2)=BB(N2)*B0BI 

00167 52 BG(N2) = Y AH 2 


+GGIN3)*G1GI 
♦CB(N3T*diBl 


+GSTiiZT4* O&i 


+&6iN3}*B10I 

♦GGW>*GlBi 


00166 

00169 

0C17C 


♦Bo1N3)*6161 
♦BB(N3)*BlbI 


♦ BB (N2)*bOgI 
♦BG1N2)*G0BI 


♦ BB(N3)*B1GI 

♦ BG(N 3)*G1BI 


YAH3 


zSftlll * 

&&1I1 —*= 


=££tU*GQGX 
6B(11 *B08i 
YAH 3 


♦ bb( 1) * BOG I 
♦ Gb(ll *0081 


0C171 

00172 

00173 


YAH* 
BBC 1) 
BG( 1) 


BG ( 1 I 
= BUll) 
* YAH* 


♦GOG I 
♦BOB J 


BBU) 
BG (1) 


*B0GI 
♦ GOBI 


51 CONTINUE 
90 CONTINUE 


{FINEST.01.33) PRINT 503 
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00177 5 03 FORMAT ( / // // // ////// /50 X* 33HI SYHBUL DA T A STARTS UN NEXT PAGE ) ) 

00178 IFtNbST • GT.35) PRINT 502 

00179 SO2 FORMAT (1H1) 

06180 ' T»klNl 506, TIMEX 

00181 500 FORMATI lx,A8/ IX, SSYMBUL NUMBtftSlOX, l5HftECUft$I0N PftOBS, 

iL_17HC0NDIT TONAL PROBS*__ 15Xjr_ 13HCUM RECURSION* 12X, *1. - CUHS, 

2 9X• SMEAN •TWEENS) 

00182 P = AlPHA*(Gb(l> ♦ GB(1)) ♦ BETA * { BG ( 1 > ♦ BB(1)> 

00183 _ Q « P _ 

00184 ft 1 * !♦ — P 

00185 51 * Q * FMUO / R1 - 

00186 PRINT 507* P, 0* Rl, SI 

00187 507 FORMAT (12X,$0*$, 8X, E16.9, 1H*, 34X, E18.12, 1H*, 9X, 

1 E10 • 4 , S*S, 9X, b10 • 4 * S*S) 

00188_ P BAR = 0 -____ 

00189 PYAR = _ 

00V9G DO 53 I-2,ft5 


00191 1S _ « I - 1 

00192 P = ALPHAMGGU) ♦ GB ( I ) ) ♦ 8 E T A *(BG(I> ♦ BB ( I ) ) 

00193 Q = P ♦ 0 

00194 _ QH = 1. - Q _ 

00195 R = P / R1 

00196 S = 0 * FMUO / QH 

00197 FIS * FLOATUS) 

00198 FIST* FIS**2 

00199 P BAR = PBAR ♦ P * FIS 

00200 PVAR = PVAR ♦ P * FIST 

00201 _ PRINT 501, IS, P, Rf 6V OH, 5 

66202 53 CONTINUE 

00203 501 FORHA?_t 8X*_fS, \H* f 21 9x, £15.8, !«♦), 9X, E18.12 , |H*, 9Xt 


1 E10.4, 1H* , 9X, E10.4, $*S> 

00204 606 CONTINUE 

00205 VAR * P VAR - PBAR*PBAR 


00206 

66207 

00208 

00209 

00210 


06211 

00212 

00213 

00214 


SVAR 


SORT < VAR) 


CMEAN = PBAR / ftl 

CAU. Tf*£ 


PRINT 305, PBAR, 
305 FORMAT (// 25X, 

1 25X 


TIMEX 


3000 CONT INUE_ 

if(IK .EQ. 
GO TO 404 
ENO 


25x, 

25X, 


VAR, SVAR, CMEAN 
SPARTIAL MEAN 
SAPPROXIMATE VARIANCE 
$APPROXlMATE ST. DEV. 


= S, 


E13•6/ 
El 3.6/ 


1, 

SPARTIAL CONDITIONAL MEAN * S, E13.6//1X, A9//I 


0) GO TO 604 
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